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Abstract 

We study the generation of gravitational radiation by sources moving in 
the de Sitter background. Exploiting the maximal symmetry and the con- 
formal flatness of de Sitter space-time we prove that the derivation of this 
gravitational radiation can be done along the same lines as in Minkowski 
space-time. A gauge is chosen in which all the physical and unphysical modes 
of the graviton are those of a minimally coupled massless scalar field in de 
Sitter space-time and a massless field in Minkowski space-time. The gravi- 
ton retarded Green's function and the Schwinger commutator function are 
computed in this gauge using Quantum Field Theory techniques. We obtain 
closed formulae for the spectral decomposition in frequencies of the linearized 
gravitational field produced by the source, in terms of a suitable spectral 
decomposition of the source energy-momentum tensor Tjp u . This spectral 
decomposition is dictated by the free (sourceless) gravitational wave modes 
in the de Sitter background. 

04.30.-w; 04.30.Db; 04.62,+v 



Typeset using REVTeX 



*Laboratoire Associe au CNRS UMR 7589. 



1 



I. INTRODUCTION 



De Sitter space-time is interesting for several reasons. To begin with, the de Sitter space- 
time is the natural framework for inflationary models of the universe, in which the universe 
experiences a period of exponential expansion driven by the energy density of the vacuum. In 
addition, it is the simplest model of space-time with a non-vanishing cosmological constant 
A; i.e. it is the simplest model of empty space-time but for the energy of the vacuum. 
It shares with Minkowski space-time the property of being maximally symmetric, but on 
the other hand it has a non vanishing curvature. Thus, it provides a tractable example of 
space-time in which the effects of curvature can be explicitly computed. 

Gravitational radiation is a subject of fundamental importance and current interest in 
both astrophysics and cosmology. Its direct detection is now a realistic challenging pos- 
sibility within the current and expected level of detector's sensitivity. Besides the weak 
field sources (i.e. binary star systems, fissioning stars, oscillating and rotating spheroids,...), 
neutron stars, collapsing supernovae, quasars and black hole collisions are sources of intense 
gravitational radiation. In cosmology, gravitational waves would arise during the several 
phase transitions undergone by the early universe, as well as from string sources, and at 
the end of inflation []]]. The linearization procedure around flat Minkowski space-time is 
clearly not applicable in such strong field problems. Perturbation techniques around curved 
backgrounds, and in the presence of sources, are needed. 

In this paper we consider the generation of gravitational waves by generic sources which 
move in the de Sitter space-time. That is, we compute the waves produced by a source 
described by an energy-momentun tensor T^ u , which is covariantly conserved with respect 
to the de Sitter metric. It turns out that the maximal symmetry of the de Sitter space-time, 
and the use of conformal coordinates to exploit its conformal flatness, allow a discussion of 
gravitational radiation in this space following analogous steps as in Minkowski space-time. 
Of course, we consider weak enough sources for the linear approximation to hold. Apart 
from this condition, the sources are totally generic. 

As one of the main results of our study, we obtain explicit closed formulae for the 
spectrum -amplitude as a function of the frequency- of gravitational waves in terms of a 
suitable spectral decomposition of the energy-momentum tensor which is producing these 
waves. These formulae should be interesting for primordial cosmology, since they relate the 
spectrum of primordial gravity waves produced during the inflationary period to the energy- 
momentum tensor of the possible sources existing during this period. This is a classical 
mechanism of production of primordial gravity waves, and it is different to the quantum 
mechanical production of gravitons due to the time variation of the de Sitter metric which 
has been studied before M, H. 

The paper is organized in a self-contained way as follows: In section II we consider 
a linearized gravitational perturbation around the 4-D de Sitter background. Thus, we 
consider the pertubative expansion of the full Einstein equations -including the source and 
the metric perturbation- in orders of the perturbation h 1 ^. Due to the maximal symmetry of 
the de Sitter space-time, the Einstein equations up to first order, can be consistently splitted 
into two equations: a zero order equation describing the de Sitter space being produced by 
the cosmological constant, and a first order equation which describes the production of 
gravity waves from a source. Thanks to the maximal symmetry, the covariant conservation 
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of the source energy-momentum tensor T^ u with respect to the de Sitter metric coincides 
with the integrability condition for the first order part of the Einstein equation. 

In section III we discuss the gauge invariance of the first order equation. Again the 
maximal symmetry of the de Sitter space-time implies the gauge invariance of the first order 
part of the Einstein tensor G^> v under infinitesimal coordinate transformations. Then, in 
order to give the first order part of the Einstein equation a definite form, we set a gauge 
fixing condition of the form D^tp^ = B^, where i\f is the trace reversed graviton field. The 
field is chosen as to eliminate the non-diagonal terms in first derivatives from G^ u . With 
this choice, we arrive to the simple decoupled system of wave equations ( |3.23| ) governing the 
linearized gravitational field ip v produced by the source Tjf> v . 

In section IV we solve the homogeneous version of eq.( |3.23 ) which describes free (source- 
less) gravitational waves. The solutions turn out to be remarkably simple as a consequence 
of our gauge choice ( 3.21|) . They amount to de Sitter minimally coupled massless scalar 
modes and Minkowski massless modes ||. We also discuss the residual gauge invariance 
allowed by our gauge choice and use it to extract the two transverse traceless physical po- 
larizations of the graviton. We conclude this section setting the conditions of validity for 
the linear approximation. 

Section V is devoted to the computation of the retarded Green's function solving the 
graviton wave equation Q3.23| ). The graviton retarded Green's function and also the Feynman 
propagator in de the Sitter space-time have been studied in previous papers 0,0, following 
a geometrical approach developed in M. In this paper we compute the retarded Green's 
function most easily using QFT techniques. By the way we also obtain a very simple 
expression for the Schwinger commutator function of a minimally coupled massless scalar 
field. 

Finally in section VI we compute the linearized gravitational field produced by a generic 
source with energy-momentum tensor T^ u . We show that the spectral decomposition of 
this gravitational field in frequencies, can be easily expressed in terms of a suitable spectral 
decomposition of the energy momentum tensor T^ u . This spectral decomposition is dic- 
tated by the form of the free gravitational wave modes. We also show that the linearized 
gravitational field produced by localized sources takes the form of free gravitational waves 
being radiated away from the source, and whose amplitudes are easily related to Tjf> v ' . 

We end up with the conclusions in section VII, followed by an appendix with the formulae 
for covariant derivatives, curvature tensors and d'Alembertians for the de Sitter space-time 
which are needed in the paper. 



II. GRAVITATIONAL PERTURBATIONS IN DE SITTER SPACE-TIME 

In this section we study the linearized gravitational perturbations around the 4-D de 
Sitter background. Thus, we start from the metric 

SV = 7m" + V* = JpT~2 + <?V) (2- 1 ) 

where 77^ = diag (—,+,+, +) is the Minkowski metric, 7^ is the de Sitter metric, and 
h^y = (H 2 r] 2 )~~ 1 <j) f j iL/ is a small perturbation, in the sense that the components of the tensor 
density (f)^ are much smaller than 1. 
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As it is less involved to work with tensors rather than with tensor densities, it is conve- 
nient to introduce the tensor field 



K = ^Kp (2-2) 
where r ) vp is the inverse de Sitter metric, and so 

/*o = -0oo , /i? = = , K h) o u (2.3) 
It is also convenient to introduce the trace reversed graviton field 

% = K~\^ h (2-4) 

where h = h\. 

A remark about the conventions for the raising and lowering of tensor indices is here in 
order. The indices of h 1 ^ (or ip'p are raised and lowered with the de Sitter metric, 7 Miy denotes 
the inverse of 7^, denotes the covariant derivative with respect to the de Sitter metric, 
and D u = r ) vp D p . For the rest of the tensors, their indices are raised and lowered with the 
full metric and its inverse g pu , which up to first order reads g pv = 7^ — h pu + 0(h 2 ). Notice 
however, that when a tensor has a perturbative expansion, the former rule applies to the 
full tensor, and not to each of the terms in the expansion. Thus if we consider, for example, 
the expansions for the full Ricci tensor up to first order in h v ^. R p = R^ u + R^ u + ■ ■ ■ 
and R pK = + Rfy + •••■, we have R v = g up R Pf)) which implies for the first order term 

Now we consider the Einstein equations 

g^r^-\^r = -mqt; (2.5) 

and deal with them in the following way: we develop the Einstein tensor G v up to first order 
in h v u 

G; = G$> + G%> + ■■■ (2.6) 
and for the second member we set 

T; = T { ^ )v + T^ v + ■■■ (2.7) 

where 

Tf>» = -AS» (2.8) 

is the energy- momentum tensor producing the de Sitter background with cosmological con- 
stant A = and Tjp v is the energy-momentum tensor of a source. This source is moving 
in the de Sitter background. So, T^ u is covariantly conserved with respect to the de Sitter 
metric, i.e. 

D v Tjp v = (2.9) 
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Thus, up to first order in the Einstein equations Q2.5Q give the two equations 

G i f v = -MTf> (2.10) 

and 

G<$ v = -%-kQT$> (2.11) 

Then, equation ( p.lOj ) describes how the background is produced by the energy- 
momentum tensor Tjf^, and ( |2.11 ) describes the production of a gravitational perturba- 
tion by the energy-momentum tensor T^ u , which is to be considered small, i.e. of first 
order. This is exactly the point of view adopted when computing gravitational radiation 
by astrophysical sources in Minkowski space-time |§, [10|, with the only difference that, for 
Minkowski space-time, T^ u vanishes. There is however a catch: we should make sure that 
eq. (|2.11|) is consistent, i.e. integrable. Now, the integrability condition for the full Einstein 
equation ( |2.5| ) takes the form 

V V T; = D v Tf> + D V T^> + - rg]" Tf> + • • • = (2.12) 

where is the covariant derivative with respect to the full metric g^. In eq. Q2.12 ), we 
have developed the metric connection up to first order as 

r; A = r^ + r^ + --- (2.13) 

with being the metric connection for the de Sitter metric given in ( |A3| ) , and its 
first order correction 

rg" = \ {d,K + D x h; - D-h*) (2.i4) 

Now the zero order part of eq. Q2.12| ), D U T^ U = 0, is obviously fulfilled, and the first 
order part yields the integrability condition 

D V T® V + r£j> T^ p - T^p T(° )u = (2.15) 

However, due to the maximally symmetric form (|2.8|) of T^ u , the terms proportional to 
I^ 1 ) identically cancel each other, and we are left with eq. ( |2.9|) as integrability condition 
for eq. ( |2.11|) . Thus, equation ( [2.1 1|) can be consistently solved for any de Sitter covariantly 
conserved source T^ u , and it describes, as we shall explain in this paper, the production of 
gravitational radiation by this source. In particular, for T^ u = 0, we have the equation 

Gj )v = (2.16) 
describing the propagation of free gravitational waves in the de Sitter space-time. 
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III. GRAVITATIONAL WAVE EQUATION AND CHOICE OF THE GAUGE 



Once we have set equation (|2.11|) , our next step in order to solve it, will be to give it a 
definite form. This entails a combination of two things: to obtain the expression of 
in terms of and to fix the gauge. For the first, we start with the expression of the Ricci 
tensor in terms of the metric connection 

R, K = d K v^ x - d x v* K + r£ A r* p - r£ K v x pX (3.1) 

and expand it up to first order 

R, K = RQ + BQ + ■■■ (3.2) 

where Rj® is given in and 



RM = - 



D X D X V + D K D^h - D X DX - D x D K h* 



(3.3) 



Then the l.h.s. of ( p.llj ) is given by 



(%> = B$> -U; Rp) = - h™R$> - I 5; RM (3.4) 

with i*M = R[ 1)x 

Next, we come to the gauge fixing. As a consequence of the covariance of the Einstein 
equations ( |2.5| ), and the maximal symmetry of the de Sitter space-time, the equation for the 
gravitational perturbation fl2.11|) is gauge invariant under the transformations 

— ► h'^ = hfut + D^ K + D K ^ (3.5) 

where £ M is an infinitesimal vector field to be considered of the same order as h v , i.e. first 
order. Let us explain how this gauge invariance arises. Under an infinitesimal coordinate 
transformation 

x » — > x '» = x fi - e(x) (3.6) 

the metric g^ u transforms as 

9 iv — * 9'^ = 9^ + C-n g^v = g^u + D^ v + Dvtn (3.7) 

where £g is the Lie derivative for the vector field and we have omitted second order 
terms. Thus, we can take the point of view in which the background metric 7^ remains 
invariant under the infinitesimal coordinate transformations (|3.6|), while the perturbation 



h^ v transforms as stated in ( |3.5| ). Then, for the tensors G^ and T" the zero order parts 
remain invariant under (|3.6|), while the first order parts transform as 



GW"— >C«" = GjJ> + £ t G£> (3.8) 

and a similar expression for T^ v ' . But now, the maximally symmetric forms G^ u oc 5^ 
and T^ v oc 6", imply C ( G^> = and = 0. Thus, is invariant under the 
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infinitesimal coordinate transformations Q3.6|) , as can also be directly checked by substitution 
of the gauge transformation (|3.5| ) in the expression (|3.4j) . Hence, eq. ( p.ll| ) is a second order 
partial differential equation for the unknown function h 1 ^, with a gauge invariant source 
T^ u , which is gauge invariant under the transformation fl3.5| ). Therefore, in order to solve 
it, we need to set a gauge fixing condition on h v (or ip"). Of course, if we were to solve 
only the homogeneous equation ( |2.16| ), the simplest procedure would be to impose the gauge 
condition [01, H. 



(3.9) 



and then use the residual gauge invariance allowed by (|3.9|) , for going the so called TT 
gauge [ItJ, by means of the additional condition tj) = 0. However, this cannot be done 
in the presence of a non- vanishing source term T^ u . The situation here is similar to the 
Maxwell equations with a source term: OA^ — d^(d v A u ) = J M . Once we have chosen the 
Lorentz gauge d v A v = 0, we cannot impose in addition A = 0, unless J — 0. It is true that 
we could set indeed the gauge condition fl3.9|), and then proceed to solve the resulting gauge 
fixed form of ( |2.11| ). However, as we shall explain below, this introduces spurious (gauge) 
complications in the mode solutions to ( |2.11| ). Thus, we shall impose the more general gauge 
fixing condition 



(3.10) 



and let the equations choose their favourite gauge fixing field instead. (Notice that 
eq. ( |3.10|) implies that the gauge fixing field B^ is to be considered of first order in the 
perturbation expansion). 

Imposing the gauge condition ( |3.10|) , the first order part of the Ricci tensor (|3.3j )can be 
written as 



1 r 



D X D X h, K + 2 /^<L A + h^Rfl + KR^ K - D K B, - D,B K 
Then from ( |3.11|) and ( |3.4j ), the first order part of the Einstein tensor is 



(3.11) 



= \ d"d x r» + Rfl; h x ° + \s; i?s 



-d ,/ b^-d^b u + 5;d x b x 



Now according to ( |A7|) , the gauge fixing condition (|3.10|) reads 



" T) T) ^ 



(3.12) 



(3.13) 



This allows us to rewrite the two terms in the second line of the tensor d'Alembertian ( |A13| ) 

as 



(3.14) 
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So, for a tensor field that satisfies the gauge fixing condition ( p,13|) , the tensor 
d'Alembertian can be recast as 



H 2 r] 



2 r 
+ - 

V 



V 



„/.!> 



„/.0 „/.0 



In addition, from eqs. flAq), (|A9|) and ( |2.4|) , we have 



2 ""-Acr/i 



and 



1 *2^ = ^U 



r/- 



(3.15) 



(3.16) 



(3.17) 



Also, using the expressions (|A4| ) and ( |A5|) , the second line in eq. ( |3.12|) takes the form 
1 



-D V B^-D,B» + 5:D X B X 



-rf a d,B a - rT d a B, + 51 ^ d a B p 



Finally, adding up the expressions ( |3.15| ), (|3.16|) , ( |3.17| ) and fl3.18| ), we obtain 



(3.18) 



W = 

M #2^2 



2 



+- 

This expression suggests a suitable choice for 5^, namely 

B» = 2Hu v % 



(3.19) 



(3.20) 



where u u = Ht]5q is a unit time-like vector in the background. With this choice, the non- 
diagonal terms in first derivatives in the first line of eq. ( |3.19| ), cancel against the first two 
terms of the third line. Moreover, with the choice ( |3.20| ), the gauge fixing condition (|3.10|) 
takes the form 



D x ift - 2H u x ^ = d x ^--^ + -6°^ = 
Then, using eqs. (|3.20|) and ( [3.21|) , we also have 



(3.21) 



r, XK d x B K = ^-[^-il; 



(3.22) 
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which produces a further cancelation of two terms in (|3.19|) . Therefore, we arrive at the 
following final gauge fixed form for the gravitational wave equations with a generic source 
term 



16?rg 

H 2 V 2 M 



(3.23) 



Notice that this wave equation -which we have obtained by imposing the gauge condition 
( p.21| )- is much simpler than the one we would have obtained imposing (|3.9| ), i.e. setting 

= in eq. fl3.19j ). The reader can also directly check that the solutions to ( |3.23| ), satisfy 
the gauge condition ( |3.9| ) provided T^ u is de Sitter covariantly conserved. 

Once we have obtained the graviton wave equation ( [3.23D , our next task will be to 
solve it for a generic T^ u . For this purpose, we need the retarded Green's function for 
the differential operator in the left hand side of eq. fl3.23|) . This Green's function will be 
computed in section IV as an appropriate superposition of solutions to the homogeneous 



version of eq. (|3.23| ). Thus, we discuss first the solutions of the homogeneous equation, 
which are interesting in their own right, since they represent the free gravitational waves in 
de Sitter space-time. 



IV. FREE GRAVITATIONAL WAVES IN DE SITTER SPACE-TIME 



In this section we set the source = 0, and solve the homogeneous version of ( 3.23 ). 

2 _ ... 2 

7] "~ 7] 







(4.1) 



We shall see that in the gauge ( |3.21|) we have a very simple basis of mode solutions to 
the free gravitational wave equations Q4.1|) . To describe the solutions of (|4.1j) it is convenient 
to introduce the tensor density 



whose components are related to those of ty v by 

Xoo = -V>o > Xoi = -tpi = , Xij = tpl = i>) 
In terms of Xfw the metric perturbation (p^ u is 

(fi^u = X»v - \ Vfiv Xaf3 

We shall also use the notation 

X = Xu + X22 + X33 , X = Xoo + X 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Then, splitting eq. (|4.1|) into its time-time, time-space, and space-space components, we 
obtain 



2 2 
D Xoo + - drjXoo 2 X = 



(4.6) 
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2 o 2 

n XOi + ~ Or,X0i o XOi 

7] T] z 



(4.7) 



□ + - d v j Xij = 



(4.8) 



From eqs. ( [A 10]) and ( [4.8D , we see that the space components Xij behave as minimally 
coupled massless scalar fields in de Sitter space-time. As it is known ||, the two physical 
TT components of a free gravitational wave in de Sitter space-time -which count among the 
Xij- behave as minimally coupled massless scalar fields. As we mentioned above, this result 
can be easily obtained in the absence of sources by sequentially imposing the two gauge 
conditions D u ip^ = and = 0. Then, one good property of the gauge condition ( |3.21| ) - 
which can be applied in the presence of sources- is that it alone is sufficient to capture the 
minimally coupled massless scalar field behaviour fl4.8|) . In addition, the two other equations 
( |4.7| ) and ( |4.6| ) are very easy to solve. In fact, eq. ( |4.7|) can be rewritten as 







Moreover, adding eq. 



□ \-Xoi 



) and the trace of eq. Q4.8|) we also have 
'1 N 



(4.9) 



□ 



Jl 



■x 







(4.10) 



Thus, the time-space components xoi and the combination x behave as free massless 
Minkowski fields rescaled by a factor rj. 

Now it is very easy to write down the general solution for x^v as a superposition of plane 

waves 



Xw (V, 



d 3 k 



fnv (77; k) exp (ikx ) + f* u (77; k) exp (-ikx 



For xoi and % we obviously have 



Xoi (V, x) 



d 3 k 



e 0i (kj 7] exp (—iuorj + ikx^j + e* Qi (kj rj exp (iurj — ikx 



(4.11) 



(4.12) 



and 



, ( , xl - / d 3 k [. , exp + nr) + r (,-) , e XP H - Hr) ] ,4,3, 

where ui = \k\. 

In addition, the ordinary differential equation for the modes fij(T]',k) can be solved in 
terms of Hankel functions of index 3/2 So, we have 



Xij (?7, x) = J d 3 k Cij (jPj yq j exp (—iurj + ikx ) + e*j (JPj yq H j exp {iujr] — ikx 

(4.14) 

And from equation ( |4.5|) and (|4.13|) 
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Xoo (77, x) 



d 3 k { 



+ 



(Jcj 17 — e (pj yq J exp (^—icor] + ikx^j 

e* (JcJ r] — e* (jPj ^77 H — ^ exp (2^77 — z/cx ) } 



(4.15) 



where we have defined e(fc) = d^-e^ (&). 

Eqs. (|4.12|) - (|4.14f ) provide the general solution for a free gravity wave in de Sitter space- 
time in the gauge ( |3.21| ). As we shall see in the next section, the simplicity of these solutions 
will allow us a direct derivation of the graviton retarded Green's function. In this respect, it 
is enlightening to point out, that if we had imposed D v i\) v = as our gauge fixing condition, 
we would have obtained a much more difficult coupled system of partial differential equations 
for Xij.v than eqs. fl4.6|)-( P~8|) . In particular, the solutions to those equations involve Bessel 
functions of index v = V33/2. 

In the rest of this section, we show how the physical degrees of freedom for the gravi- 
tational waves are contained in expressions (|4.12| ), (|4.13|) , and Q4.14 ). Due to the linearity, 
it is enough to consider just one mode with wave vector k. To begin with, the gauge fix- 
ing condition ( |3.21|) reduces the ten polarizations that we have in Xnv down to six. More 
precisely: for the field Xnv, the gauge fixing condition ( |3.21| ) reads 



d v Xoo + diXoi + ~ X = 
77 



d v Xoi + djXji + - Xoi = 
77 



(4.16) 
(4.17) 



which translated to the modes fn V (Tj)k) exp (ikx^ gives the following constraints among 
polarizations 



and 



rueoi + e - e = 



e 0i + rijCji = 



(4.18) 
(4.19) 



where n, = ki/u) is the unit wave vector. 

Next, the residual gauge invariance in the gauge fixing ( |3.21| ) allows to reduce the six 
independent polarizations to the two physical polarizations. Under the gauge transfor- 
mation (|3.5| ), the trace reversed graviton field ip v transforms as 



(4.20) 



Then, the invariance of the gauge condition Q3.21 ) under the transformation ( |4.20| ) requires 
^ to be a solution of the equation 



2H(u„D x t-uxDJ A -u x D% 







From eqs. 



, and (|A12| ) it reduces to 



(4.21) 



(4.22) 
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which yields the simple decoupled equations 



and 



(4.23) 



(4.24) 



The solutions to eqs.( }1.23| ) and ( 4.24j ) decompose in modes in the form 

e (,, f) = / dH, [te (S) , exp (-**, + ikx) - i el (k) „ exp - ikx) } (4.25) 



and 

e (t?, f ) = f d 3 k 



i e. 



(jcj (j) ^ exp (—iur] + ikx ) — z e* (jtj ^77 H — ) oxp ( /a.'// — ik.r 



uJ 



On the other hand, from the transformation law ( 4.20|) we find 



(4.26) 



(4.27) 



(4.28) 



6 i x = ^o v e--e 

7] 



(4.29) 



Then, from eqs. ( |4.27| )-( |4.29| ), the residual gauge transformations given by ( (4.25|) -( i.26D , 
amount to the following transformation laws for polarizations 



e' = e + AuoEr, 



eoi 



'0i 



&0i ~ kiSQ — USi 



dj — ► = eij + kiEj + kjEi + (we - kej 



(4.30) 



Thus, it is possible to eliminate the polarizations and e by suitably choosing the 
parameters e M . Therefore, we can set the supplementary gauge conditions 



Xoi = X = 



(4.31) 



Of course, the gauge transformations (|4.30|) are compatible with the constraints (|4.18|) and 
( |4.19| ), and by replacing the supplementary conditions eo« = e = in eqs. (|4.18|) and (|4.19|) 
we find 



(4.32) 
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This leave us with the traceless transverse (TT) graviton physical polarizations. Notice that 
since e = 0, eq. (|4.15|) implies Xoo — 0. Therefore, (p^ v coincides with x^u, and the TT 



conditions are satisfied in this gauge by the metric perturbation itself. 

We conclude this section with a comment about the validity of the linear approximation 
we apply. Let us consider just one physical mode for 0^, with amplitude A and wave vector 
k. 

A^r] ^ exp ^—iur] + ikx ) + A* ^77 H J exp (itur/ — ikx ) (4.33) 

As stated at the beginning of section II, the validity of the linear approximation requires 
|</V«/| "C 1. Thus, according to the functional form of the mode ( }L33|) , the amplitude A must 
satisfy the conditions 

|4|<A , \A\<^u (4.34) 

\m 

that is, the linear approximation holds for high frequencies and small conformal time. Let us 
remind that the cosmic region for the de Sitter space-time which describes an exponentially 
expanding space-time, corresponds to — 1/H < 77 < 0. Thus, the linear approximation holds 
throughout this region provided that 

A < max {if, uj} (4.35) 



V. RETARDED GREEN'S FUNCTION FOR THE MINIMALLY COUPLED 
MASSLESS SCALAR FIELD AND THE GRAVITON 



We come back now to the inhomogeneous equation fl3.23j ). Splitting (|3.23| ) we have 



□ Xoo + - «9,Xoo - ^ x = -16 TT 6? T$ (5.1) 

□ Xo* + - d vX oi - 4 Xoi = -16 TT g T$ (5.2) 
(□ +-d v )xi j = -16nGTjp (5.3) 



77 

As in previous section, eq. ( |5.2j ) can be rewritten as 



□(±xb.)=-^i£> (5.4) 

V 77 I 77 



and adding up eq. (|5.1j ) with the trace of ( |5.3|) we also have 



□ -X = -T {1) (5.5) 

\V J V 
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Here 

j. (i) = T w + T (i) + T g) ; ^ = T (%) + y(i) (5 _ 6) 

Thus, the retarded solution to eqs. ( |5.4p and ( |5.5| ) can be obtained using the well known 
Minkowski massless retarded Green's function 

Gf ] (x, x') = 47r| g_ g| * (V ~ rf - \x - f | ) (5.7) 

which has support on the past light cone and satisfies 

□ Gf ] (x, x') = -5 (4) (a; - x') (5.8) 



On the other hand, in order to solve ( |5.3| ) we need the retarded Green's function for the 
scalar d'Alembertian in de Sitter space-time, i.e. we need the retarded solution to 

D\ G R (x, x') = -#y 5 (4) (x - x') (5.9) 



Here D\ is the scalar D'alambertian in de Sitter space-time given in ( |A10|) . Thus, eq. (" |5.Q| 
takes the form 



+ G R (x,x') = -H 2 v 2 5^( 



x-x') (5.10) 

This retarded Green's function G R (x, x') could be obtained as in 0, using the geometric 
techniques developed in ||. However, although the retarded Green's function is a purely 
classical object, it can be most easily obtained using QFT techniques. Thus, we consider a 
scalar field 

4>(r],x) = J d 3 k u^(r/, x) a(k) + Ug(r], x) a)(k) (5.11) 
obeying the homogeneous equation 

^□ + ~dj0 = O (5.12) 

i.e. 4>(rj, x) is a minimally coupled scalar field in de Sitter space-time. Then, from the 
previous section, the modes (r], x) take the form 



(2tt) 3 / 2 V ' uo 
These modes satisfy the normalization 

(u%, u % ) =6{k-k') (5.14) 



with respect to the scalar product 

i 



>2 , 9\, 



H 2 V 2 



d 3 x 4>l(r], x) d v 0i (r/, x) (5.15) 
14 



Now, if we canonically quantize the scalar field according to 

~a(k) , a f {k')]= 5{k - k') (5.16) 
its retarded Green's function is given by 

G R {x,x') = -9(rj- r)')G{x,x') (5.17) 

Here 9{rf) is the Heaviside step function, and G(x, x') is the Schwinger function for the field 

(j)(r],x) 



G(x,x') = -i (0| [<f>(x), (j)(x')\ |0) 



(5.18) 



Then, by inserting the mode expansion ( 5.1 If) in ( 5.181) , a straightforward computation 
yields 



G(x,x') 



H 2 



57T Z \X — X 



- (777/ h-i(rj- 7]') I 2 + h) 



where 



f°° du 

In = / ( ex P( m 2/-) - exp{iuy+)) ; n = 1, 2, 3 

J — oo 11 



and 



y-b = r) — r\ ± |af — a?] 



(5.19) 



(5.20) 



(5.21) 



The integral l\ gives 

A = 2tt (<%_)-<%+)) (5.22) 

The integral I 2 is clearly convergent, and its evaluation using residues theorem yields 

h = 2m (0(y_) - %+)) (5.23) 

Notice that the integral ^3 is apparently logarithmically divergent in the infrared -the 
integration variable being u = \k\- since the integrand behaves for small u as 



1 i 
— ( exp(iuy-) — exp(iuy + ) ) = 2 \x — x'\ + • • 



ii- 



ii 



(5.24) 



However, this divergence does not really exist because I3 can be rewritten as 

h = I' 3 + h (5.25) 



where 



/OO fj^y^ 
— ( exp(iuy-) — exp(ra?/ + ) + 2% \x — x\ sinii 
-00 u 



(5.26) 



and 
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I3 = — 2i \x — x I / — smu (5-27) 



Then I' z is convergent, and ^3 vanish -in the sense of Cauchy principal value- because the 
integrand is an odd function. Therefore, I3 = I' 3 , and its evaluation using residues theorem, 
yields 

J 3 = 7T (\y + \-\y_\) (5.28) 

Thus, the Schwinger function G(x,x') is perfectly finite in the infrared limit, and replac- 
ing expressions (|5.22|) , (|5.23|) , and (|5.28|) in eq. (|5.19|) , it takes the form 



G(x,x') = - - ^ + \rjrj' (<%_) - 5{y+)) + (77 - rf) (%_) - %+)) 



Air \x — x' 



+ \(\v+\-\V-\)] (5-29) 



Finally, from eqs. ( |5.17| ) and ( |5.29| ) we obtain the following simple expression for the de 



Sitter scalar retarded Green's function 

G R (x, x') = - H l r]V ^ l 6(ri-rf-\x-#\) + ^0(ri-rf-\x-&\) (5.30) 
Ait \x — x'\ Air 

Notice that to obtain this expression, some cancellations have been produced in multi- 
plying eq. ( |5.29p by 8(r] — 77'). This is due to the following identities among distributions 



6(r] - rf) 


Kv+) 


= 


0(77 - 4) 


Kv-) 


= Kv-) 


0(V - rf) 


o(y+) 


= 0{V ~ rf) 


0(V - rf) 


0(y-) 


= %-) 


9(t]-t]') ( \y+\ - 


\v-\) 


= 2(77-77^(77 



77 , )-2 2 /_%_) (5.31) 

In particular, the advanced variable y + has disappeared from the arguments of the dis- 
tributions entering ( |5.30| ). Thus, the retarded Green's function ( |5.30| ) has support on the 
past light cone and its interior, as it should be according to the causality of the classical 
theory. Notice also that the main difference with the Minkowski retarded Green's function 
is the term proportional to 9{yJ). This term tell us that -although the free gravitational 
waves propagate at the speed of light- in the production of gravity waves from sources, there 
is information about these sources, which propagates at a lower speed. 

Expressions (|5.7| ) and (|5.30|) completely solve the problem of the retarded graviton prop- 
agator in the de Sitter space-time. With this propagator at hand, we are going to discuss 
in the next section the production of gravity waves by a generic source. Notice that from 
now on we shall supress the upper label (1) from the energy- momentum tensor of the 
source, in order to alleviate the notation. 
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VI. GRAVITATIONAL WAVES PRODUCED BY SOURCES IN DE SITTER 

SPACE-TIME 



We proceed to obtain the gravitational field produced -in the linear approximation- by 
a generic source with energy- momentum tensor T^ v . From eqs. (|5.7| ), ( |5.8|) and (|5.4j) the 
time-space components for this field are 

Xoifaz) =r]\Q-ng J - X ^f V G^\x,x')T 0i (x') 

= 4GV [rS^ T Ql ( V -\x-x'\,x') (6.1) 

J \X — X \ 7] — \X — X \ 

Similarly, for the tilde component we have 

X(v,x) =4:Gr] f J* 1 f{r}-\x-x!\,x*) (6.2) 



On the other hand, for the space-space components, eqs. (p.3|), ( |5.10| ), and ( |5.30| ) give 



Xij (V>x) = 16tt£ J d ^ 2 ^ G R (x,x')T i:j (x') 



V 

Then, using the identity 



d 3 x' r] 
T — x'\ rj 

4 g J d 3 x' drj \ 6 (rj - rj - \x - x'\ ) (r)', x') (6.3) 



4 Q [ — l S(ij — ri , — \x — x'\) TiAri' } x' 

J \X — X'\ vj 



— T=i S (t] — rj' — \x — x 1 ] ) -\ 6 (rj — rj — \x — x \ ) 

rj \x — x'\ ' rj' 1 



1 (I 

— — 5 (rj — rj — \x — x'\ ) — dw — - 6 (rj — rj 
— x'\ \rf 



\x-x\)\ (6.4) 



and doing a partial integration, eq. (|6.3|) can be recast as 

dh 



Xij (V,%) = 4 £ J ~ \x-x\ , f ) 



p-\x-x'\ drf 



+ 4G d'x' -L d^T^r,' x') (6.5) 

o T] 



Eqs. (|6.1|), (|6.2| ), and ( |6.5D give the linearized gravitational field Xta> produced by a source 
with energy-momentum tensor T^ v in de Sitter space-time. With this formulae at hand, we 
undertake now a double task: to do the spectral decomposition of the gravitational field Xnu 
in frequencies, relating it to the spectral decomposition of the source T^, and to show that 
this field takes the form of radiated free gravitational waves in de Sitter space-time, when 
we go to the "wave zone" far away from the sources. For the first, we need an appropriate 
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decomposition of the energy momentum tensor in modes of frequency uj. For the space-space 
components, it is suitable to decompose the energy-momentum tensor in the form 



/•oo 

T i:j (Tj,x) = du 



exp(-iiuri) rj 



UJ 



%j(u, x) + exp(iur]) (r) + — ) 7$(u, x 



UJ 



J du exp(—iuri) (w — — ^ x) 



(6.6) 



tational wave. Notice that we have defined T^—co, x) = T* u (u, x) as usual. Now, since the 



i.e. we decompose Tjj according to the modes entering the components Xij °f a f ree gravi- 
tational wave. Notice that we hav 

modes exp(— iur]) (rj — satisfy 

i d v exp(—iujr]) (rj ) = ujrjexp(— iurf) (6.7) 



UJ 



the integral transform ( |6.6| ) can be easily inverted in the form 



%j(u,x) 



2n uj 



drj exp(iujr]) - cL 7^(77, x) 



(6.8) 



For the other components of T^, since the corresponding modes of a free gravitational 
wave take the form rj exp(iujrj) , we do the spectral decompositions 



and 



whose inverses are 



and 



/oo 
duj exp(-iujr)) %j(uj,x) 
-oo 

_ roo 

T(r],x) —r\\ du exp(-iurj) T(oj,x) 



T 0i (uj,x) 



T(uj, x) 



(6.9) 
(6.10) 

(6.11) 
(6.12) 

(6.13) 



Then, for the space-space components of Xij, replacing the spectral decomposition (|6.6|) for 
Tij in the expression ( |6.5| ), we have 



drj 

2n J-oo rj 
1 A 00 drj 



exp(iur]) T 0i (uj,x) 



exp(iujfj) T(u, x) 



2n J-oo Tj 

Now, let us decompose the gravitational field Xnv in modes of frequency uj 



+4£ / d 6 x 



d 3 x' 



'3^ 



exp(— iujr/) exp(iuj\x — af\) [rj — \x — x'\ ) %j(u), 0?) + c. c. 



r\— \x— x'\ 



4 UJ %j(uj, x') 



Tj — \x — x'\ 



drj 1 exp(—iujrj') 



c. c. 



(6.14) 
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and doing the usual shift uj — > uj + ie, to handle the lower limit of the integral over 7/ in 
the second line of (|6.14 ), we find 

Xii\vi%) = ^G \V ) exp(— iujr]) \ —; — — 7 exp(za;|x — x'l) %a(oj, x 1 ) + c. c. (6.15) 

J \ uj J J \x — x'\ 

Similarly, replacing the spectral decompositions (|6.9|) and ( |6.10|) for Toj and T in the 
expressions ( p.ip and ( |6.2| ) for xoi an d X we have 

Xoi\vi %) = 4 G V exp(— iujr/) [ _f X exp(iu\x — x'\) T 0i (u;, x') + c. c. (6.16) 

J \x — x'\ 

and 

/d x' 
—; — = 7 exp(iu\x — x'l) T(u, x') + c. c. (6-17) 
\x — x I 

The expressions ( |6.15[ - (|6~i7|) are one of the main results of this paper. Just as in 
Minkowski space-time, these formulae relate the u frequency component of the gravita- 
tional field produced by the source T^, with the uj frequency component of T^ u itself. On 
the other hand, the r^-time dependent factors in front of the integrals in these formulae, 
exactly coincide with the 77-time dependent factors for a free gravitational wave of frequency 
uj. This is due to the appropriate decomposition of the energy momentum tensor which 
has been done in eqs. (|6.6|), (|6.9D, and ( |6.10| ). This decomposition is of course different 
of the plain Fourier transformation which is done in Minkowski space-time, since the form 
of the spectral transform is dictated by the form of the free gavitational waves themselves. 
Moreover, we can choose a source localized in a finite spatial region, and consider the gravi- 
tational field Xnv produced by this source in the "wave zone". To do it, we go far away from 
the source and take the region of points x that satisfy 

|x| > max{a , uj 2 a] (6.18) 

where a is the size of the source. In this limit we have as usual 
11 

exp(iuj\x — x'l) ~ exp(ikx — ikx) (6.19) 



\x — x'\ r 



with r = \x\, and k = ujx/\x\. Thus, the expressions ( 6.15| ), ( |6.16 ), and ( 6.1 7|) give in this 
limit 

ci 4 G ( i \ -* -* 

Xij (v, %) ~ — [V J exp(—iujr] + ikx) Ay(k) (6.20) 



4 G ~* ~* 

Xm\v,x) V exp(-iujr] + ikx) A 0i (k) (6.21) 

r 



4 Q -> -> 

X^ (v, x) ~ 7] exp(—iujf] + ikx) A(k) (6.22) 

r 
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which correspond to free gravitational spherical waves in de Sitter space-time, being radiated 
away from the source, and whose amplitudes are given in terms of the energy-momentum 
tensor of the source by 



Anv(k) = / d 3 x" exp(-ikx) 



(6.23) 



As a final remark, it is interesting to write down the spectral decomposition in frequencies 
of the covariant conservation law for the source energy momentum tensor T^ u . Using the 
expression (|A4|) for the de Sitter covariant derivative, the conservation equation D V T IIV = 
splits into the two equations 



d v T 00 + djT 0j + = 



(6.24) 



d„Toi + djTji H — Toi — 

Tj 



(6.25) 



Notice that these equations formally coincide with the gauge fixing conditions for Xpv (|4.16|) 
and ( [4.17 ). Nevertheless, the geometrical meaning is not the same for both, because T^ v is a 
tensor while Xiu> is a tensor density. Now, using the spectral decompositions in frequencies 
( p.6|) , ( |6.9|) and (|6.10 ), the conservation laws for the uj frequency component of the energy- 
momentum tensor takes the form 



% uj (T(uj, x) — 7li(u>, xu + djT j(uj, x) = 
i uj T Qi (uj, x) + djTji(u, x) = 



(6.26) 
(6.27) 



Notice that eq. ( p.27| ) has exactly the same form as in Minkowski space-time, while eq. 
( |6.26| ) is different because T(uj,x) — Th(uj,x) ^ % (u,x), (in fact there is no suitable way 
to define % (u,x)). Finally, using the conservation laws (|6.26| ) and (|6.27| ), one can check 



that the uj frequency component x%) °f the gravitational field x^u, given by ( |6.15|) , (|6.16|) 
and (|6.17|), fulfil the gauge fixing conditions (f4.16|) and (|4.17|). 



VII. CONCLUSIONS 

We have shown that the production of gravitational radiation from sources moving in the 
4-D de Sitter background can be studied along the same lines as for Minkowski space-time. 
The maximal symmetry and the conformal flatness of the de Sitter space-time are found 
to be two key ingredients in order to achieve this goal. In addition, we have shown that - 
although the general equations for linear gravitational perturbations are rather cumbersome- 
choosing the gauge ( |3.21| ), the equations for all (physical and unphysical) polarizations of the 
graviton decouple, and amount to the equations for a de Sitter minimally coupled massless 
scalar field and a Minkowski massless field. In this respect, it is worth remarking that the 
minimally coupled massless scalar field behaviour can be easily obtained, in the case of the 
physical polarizations, by imposing the traditional synchronous transverse traceless gauge 
conditions D u tp^ = u u ip^ = 0. However, these two conditions cannot be simultaneously 
imposed in the presence of a source, and if one imposes only the condition D u tp^ = 
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instead of (|3.21|) , one is led to a much more difficult coupled graviton wave equation than 
( p.23j ), whose solutions contain spurious complications. The same thing happens for the 
residual gauge invariance allowed by the gauge condition D v ip u ^ = 0. While the infinitesimal 
coordinate transformations preserving the gauge condition ( p 21| ) are given by vector fields 
whose equations for the time and space components decouple, and give very simple mode 
solutions; the vector fields corresponding to the residual gauge invariance allowed by D u ip^ = 
0, satisfy a much more difficult coupled system of partial differential equations. 

The main new results of this paper are given in section VI. In that section we have 
shown that decomposing the energy- momentum tensor of a given generic source T^(rj,x) 
in frequencies, by using a spectral transform dictated by the modes of the free gravita- 
tional waves in the curved background, we have very simple closed formulae relating the 
uj frequency component of the linearized gravitational field produded by the source with 
the transform T^uj^) of the energy-momentum tensor. We also show that for localized 
sources, the produced gravitational field takes the form of free gravitational waves in de Sit- 
ter space-time being radiated away from the source. Thus, the generation of gravitational 
radiation by sources in de Sitter space-time ressembles very closely to the same process 
in Minkowski space-time, the main difference being in the form of the energy-momentum 
frequency transform, which enter the formulae for the amplitudes of the radiated waves. 

As a previous step we have also shown in section V, how the graviton retarded Green's 
function in the Sitter space-time -needed to solve the graviton wave equation- can be easily 
obtained using QFT techniques. The most prominent feature of this retarded Green's func- 
tion is that in addition to a delta function term in the retarded time, it also contains a term 
proportional to the Heaviside step function of the retarded time. This second term shows 
that the information about the sources in de Sitter space-time propagates not only at the 
speed of light but also at a lower speed. 

In our opinion, it would be very interesting to apply the general formulae (|6.15|) - (|6.17|) 
that we have derived in this paper, to sources that could exist during the inflationary period 
of the universe. As a first example we have work in progress [12 concerning string sources, 
whose equations of motion in de Sitter space-time have been solved in the case of a ring 
ansatz IBI. 
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APPENDIX A: APPENDIX 



In this appendix we collect a number of formulae for covariant derivatives, curvature 
tensors and d'Alembertians for de Sitter space-time, which are used in this paper. The four 
dimensional de Sitter metric in conformal coordinates x^ = (r), x) reads 

ds 2 = 7,^ dx^dx v = — {-drj 2 + dx 2 ) (Al) 

i.e. 

1 

~~ H 2 T] 2 ^ 

where T]^ = diag(— , +, +, +) is the Minkowski metric, and H is the Hubble constant. 
Then the metric connection, can be written in these coordinates as 

MV _ 1 fr ™ , rf) ™ 



r K F = --\$% + %% + %v,*) (A3) 

Thus, the covariant derivatives for covariant and contravariant vectors read 

Dfj,V K = d^V K + l - (b\ V K + 5° K V P + Vim V ) (A4) 



DpV = dJT - l - (5; V° + 81 V + 51 V x ) (A5) 
In particular we have 

D X V X = d x V x - -V° (A6) 
77 

and 

D^ = d^-Ul + Ul^ (A7) 

Since the de Sitter space-time is maximally symmetric, the Riemann and Ricci tensors 
take the form 

R { '1x p = H 2 { 1 , p1kX - 1 , x1kp ) (A8) 

and 

Rf K = -3H 2 (A9) 

In addition to these tensors, we need the scalar, vector and tensor d'Alembertians, which 
can be computed using the expression (|A3|) for the metric connection. For a scalar field <p, 
we have the scalar d'Alembertian 
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H 2 r]'- 



D X D X , 



(A10) 



where □ is the Minkowski d'Alembertian 



For a vector field the vector d'Alembertian is 

! D X D\ ^ = □ ^ + - [ d^ + <? V afS d a & 



H 2 T] 



+ 



3^ + 2^£o 



Finally, for a tensor field ip^, a long but straightforward computation yields 
D X D X ^ = □ % + 2 ~ [ d v r» + - if* 

2 r 



+ - 

rj 
2 

+ - 

7] z 



% + 25" r +2%W - <W - « 



(All) 



(A12) 



(A13) 
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